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ABSTRACT 

In a search for lensed infrared galaxies, ISOCAM images have been obtained toward 
the rich clusters Abell 2218 and Abell 2219 at 15 fim. Nine galaxies (four in Abeh 2218 
and five in Abell 2219) were detected with flux levels in the range 530-1100 fj,Jy. Three 
of the galaxies detected in Abell 2218 have previously known redshifts; of these one is 
a foreground galaxy and the other two are lensed background galaxies at z = 0.474 and 
z = 1.032. One of the objects detected in the field of Abell 2219 is a faint, optically 
red, extreme infrared-dominated galaxy with a probable redshift of 1.048. 

Subject headings: Galaxies: Clusters 



1. Introduction 



Clusters of galaxies can serve as windows on the distant universe by bringing faint objects 
above detection threshholds via gravitational lensing. Many arcs representing magnified ordinary 
galaxies at moderate to high redshifts are seen in rich galaxy clusters, and these have been 
well-studied over the past several years. Indeed some of the most distant known galaxies have 
been discovered in this way at optical wavelengths (see Pranx et al 1997 for a lensed galaxy at 
z = 4.92). Clusters are known to be equally useful as cosmic magnifiers in other wavebands as 
well. 
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In the local universe the dominant populations of luminous objects are infrared galaxies and 
quasars. The ultraluminous infrared galaxies (ULIRGs), with bolometric luminosities ^ 10^^ Lq, 
appear to be powered by intense starbursts and obscured quasars (e.g., Sanders and Mirabel 
1996). It remains an open question how many such objects exist among galaxy populations at high 
redshift, since they are obscured in the optical region and generate most of their luminosity in the 
mid- to far-infrared where sensitivities of space-borne telescopes such as IRAS have been marginal 
for their detection above redshifts of a few tenths. Indeed only one such high-z galaxy is known in 
the IRAS database, the infrared galaxy/hidden quasar F10214+4724, at z = 2.3 (Rowan-Robinson 
et al 1991). Two high-redshift IRAS quasars have also been identified: the Cloverleaf at z = 2.6, 
and APM 08279+5255 at z = 3.9, both broad absorption line quasars (Barvainis et al 1995; Irwin 
et al 1998). Boosting by gravitational lensing was required for the detection of F10214-I-4724 
(Broadhurst and Lehar 1995), the Cloverleaf is a quad lens, and APM 08279-1-5255 appears to be a 
close optical double and very likely lensed as well. New Icnscd iand unlensed infrared galaxies are 
currently being found via their submillimcter emission by SCUBA on the JCMT (Small, Ivison, &; 
Blain 1997; Ivison et al 1998; Hughes et al 1998; Barger et al 1998). 

ESA's Infrared Space Observatory (ISO) offered a new opportunity to mount systematic 
searches for other high-z infrared-dominated objects. The approach we adopted with ISO, which 
we report here, was to use the mid-infrared camera ISOCAM to image at 15 //m two very rich 
lensing clusters of galaxies, Abell 2218 and Abell 2219. The strategy was to enhance detectability 
of distant infrared galaxies over random fields by taking advantage of the cluster lensing boost. 
Other groups have carried out similar programs: Altieri et al (1998a) for Abell 2218; Altieri et al 
(1998b) and Lemonon et al (1998) for Abell 2390; and Metcalfe et al (1998) for Abell 370. 

The two target clusters were chosen for their richness and for previous optical indications of 
lensing, along with purely observational considerations such as visibility to ISO and high ecliptic 
latitude. Both are at moderate redshift: z = 0.176 for Abell 2218, and z = 0.225 for Abell 2219. 
In this experiment several probable cluster galaxies were detected at 15 jim in the two clusters. 
Three background galaxies, dominated by their mid-infrared emission, were also detected. Two 
are previously known lensed galaxies behind Abell 2218 at redshifts z = 0.474 and z = 1.034, and 
the third is a faint red object in the field of Abell 2219 with a probable redshift of z = 1.048. 

2. Observations and Data Analysis 

Each cluster was observed for a total of 1.2 hours using the LW3 filter of ISOCAM covering 
the range 12 — 18 /im. The 32x32 detector array was configured with a plate scale of 6" per pixel. 
Data were taken in "micro-scan" mode, using 3x3 rasters with 14" = 2.33 pixel steps to maximize 
the area covered and the different pixel sampling of the same sky area. Three separate such 
rasters, with slightly shifted centers, were taken of each cluster, to achieve yet more cross-sampling 
between camera pixels and sky pixels. The diffraction-limited beamsize was 6.3" FWHM, but with 
a 6" pixel size and multiple sampling, we estimate that the point source width should be 8 — 9". 
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There are no strong sources in the field with which to measure an accurate point spread function, 
but the program sources that were detected with reasonable SNR yield a FWHM of 9" when fitted 
with a circular Gaussian approximation to the PSF, or about the expected value. 

For Abell 2218 the first raster was rendered unusable by detector instabilities. Only the 
second and third rasters were summed, for a total integration time of 0.8 hr. For Abcll 2219 all 
three rasters were acceptable. Images were dark-subtracted, flat-fielded, deglitched, corrected for 
transient response, and combined using standard CIA routines at IPAC in Pasadena. 

Source fluxes were derived by fitting with a circular two-dimensional Gaussian function and 
zero level using AIPS, with fixed FWHM of 9" (see above). The RMS noise levels, obtained by 
fitting with the same Gaussian at random points in the fields, are ^ 200 /xJy for Abell 2218 and 
110 /LtJy for Abell 2219. These noise levels are well above the theoretical expectation, for reasons 
that axe not understood by us at present; the cause is not a scaling error, since our source fluxes 
are consistent with those of Alieri et al (1998a) (see §3.1). 

Absolute coordinates on the original ISO rasters have an uncertainty of ~ 6", and small shifts 
were required for both ISO images to align the ISO sources with galaxies in optical frames. No 
rotations were required. After the shifts all of the ISO sources had clear optical identifications. 

One of the detected sources in the Abell 2219 image (designated A2219#5, see below) is 
optically faint and very red, making it a candidate for the sort of distant infrared-dominated 
galaxy we were searching for. A spectrum obtained for us by T. Broadhiirst and B. Frye using the 
Keck I telescope shows an emission line at A7634A, to which we assign a probable identification 
of [on] AA3726, 3729A at z = 1.048. This identification is supported by the continuum shape and 
the lack of other strong lines in the 6300 — 9500A passband observed (see, e.g., starburst galaxy 
templates of Kinney et al 1996). A detection of Ha at 1.34 /v,m would confirm the redshift. A 
possible alternative identification would be Ha at z = 0.163, but we feel this is less likely than 
[OH] at z = 1.048 for the reasons given. 

3. Results and Discussion 

The 15 /nm images are shown in Figures 1 and 2. The coordinates shown have been adjusted 
after comparison between the sources in the ISO image and galaxies in I-band optical images, 
and should be accurate to 1 — 2". In Abell 2218 there are four sources with detections for which 
we can be reasonably confident, and in Abell 2219 we identify five sources. Only detections that 
were evident on more than one frame are considered bona fide sources here. All of the accepted 
sources are closely coincident with visible- light sources. Contour overlays of the 15 /im images on 
/-band images from the Palomar 5m telescope are shown in Figures 3 and 4 (5m images courtesy 
I. Small). An overlay of the 15 fim image of Abell 2218 on an HST i?-band image is shown in 
Figure 5. Object identifications, with 15 /xm and /-band fluxes, are given in Table 1. All sources 
are consistent with being point-like in the ISO images. 
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3.1. Detected galaxies: General discussion 

Of the nine objects detected at 15 ^m, three have redshifts available in the hterature. These 
are A2218#395, #317, and #289 (using the galaxy numbering system of Le Borgne, Pello, & 
Sanahuja 1992; Abell 2219 has no previous numbering system and the ISO sources are numbered 
here according to RA order). The latter two (see §3.2) arc background galaxies at z = 0.474 
and z = 1.034 respectively (Ebbels et al 1997), with #289 clearly being distorted by the cluster 
potential. Abeh 2218 has a mean redshift of 0.176 (Le Borgne et al 1992). A2218#395, at 
z = 0.103, is a foreground galaxy and not part of the cluster. As discussed above, we assign a 
tentative redshift of z = 1.048 for A2219#5. 

The other detected galaxies in Abcll 2218 and Abell 2219 are of unknown redshift and precise 
Hubble type, although the majority appear to be spirals or irregulars. None of those galaxies 
look particularly distorted as if by lensing, and all except A2219#5 have optical colors typical 
of ordinary spiral or elliptical galaxies. However, they are distinguished by their strong 15 //m 
emission, which sets them apart from the hundreds of other galaxies in the ISOCAM fields. 

The apparent mid-infrared luminosities (over the Au/u ~ 0.25 of the filter passband) of the 
four galaxies with known redshifts, A2218#395, #317, and #289, and A2219#5, are respectively 
Luf^m = 3.9 X IO^Lq, Lio^m = 4.9 X IO^Lq, Lr^m = 3.0 x W^^Lq, and Lr^^ = 1.9 x W^^Lq 
(assuming z = 1.048 for A2219#5), where the subscripts on the luminosities represent the rest 
wavelengths observed. The frequency width of the 15 /um filter has been taken to be 5.0 x 10^^ Hz, 
and the cosmological parameters used are Ho = 65 km Mpc~^ and = O-l- For A2218#289 
and #317, and A2219#5, the calculated luminosities represent upper limits, since the fluxes are 
likely to be magnified by lensing. The other galaxies, if assumed to be at the redshifts of their 
clusters, range in luminosity between 5.8 x IO^Lq and 2.4 x IO^Lq. For comparison, of the 57 
Virgo cluster galaxies of all Hubble types detected by Boselli et al (1998) at 15 fim, none have 
above IO^Lq, and most lie in the range 10^ — IO^Lq. The galaxies detected here have 
much higher apparent luminosities in the mid-infrared than average cluster galaxies. 

Boselli et al (1998) find that for spirals and irregulars the 15 /xm flux tends to be dominated 
by dust re-emission of primarily UV stellar light, whereas for ellipticals it is the (very weak) direct 
light of the Rayleigh-Jeans tail of the old population of red stars. However, for spirals the relation 
between star formation and 15 /um flux is not simple. It appears that galaxies with moderate 
activity have the highest 15 //m to UV flux ratios whereas for active star-forming galaxies this 
ratio is lower (Boselli et al 1997). The galaxies with clear spiral morphology in our sample are 
A2218#395 and #317, and A2219#l and #3. A2218#289 and #275 appear irregular. The 
others, all in Abell 2219 (galaxies #2, #4, and #5), are difficult to classify because of lack of 
angular resolution (all are smaller than 3"). However, the very small number of detections of 
E/SO galaxies in Virgo relative to later types by Boselli et al (1998) suggests that all or almost 
all of the detections here arc likely to be from spirals or irregulars. This is supported by the large 
15 /nm to I-band flux ratios, which are typical of spirals rather than ellipticals (see Table 1). 
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Altieri et al (1998a) recently reported ISOCAM imaging of Abell 2218 at 5, 7, 10, and 15 //m. 
Their field was roughly one-quarter the size of ours, because their pixel field of view was 3" 
compared with our 6". In the inner regions of the cluster they detected the central cD galaxy and 
five others: #395, #317, #323, #373, and #275. We detected #395, #317, #275, and possibly 
#373 (which we do not claim as a firm detection, although it does show up with ~ 400 fiJy in two 
positive contours in Figures 4 & 5). A2218#323 is weak in Figure 2 of Altieri et al, and appears 
to be below our detection threshhold, as does the cD galaxy. For #395 and #317 the fluxes given 
by Altieri et al (1998a) are consistent with ours. Our approximate flux of ~ 400 /uJy for #373 is 
also consistent with the value found by Altieri et al. 



3.2. Notes on individual galaxies 

A2218#289, at z = 1.034, consists of a complex of bright knots and diffuse emission (see 
Figure 5). Lensing distortion is substantial at the northeastern end, where the galaxy is stretched 
across the halo of cluster galaxy #244; because the object appears so luminous it is probably 
highly magnified, according to Kneib et al (1996). 

A2218#317 is a background spiral galaxy at z = 0.474. It is very probably magnified by the 
cluster, since its shape and color led to a redshift prediction by Kneib et al (1996), using cluster 
inversion techniques, of 0.2 < ^lensing < 0.4. Given that such redshift predictions are statistical in 
nature, the spectroscopically measured redshift of z = 0.474 can be considered to be consistent 
with the lensing hypothesis. 

A2219#5 is optically the faintest object among those detected at 15 /um, by more than an 
order of magnitude. It is quite red, with an optical color B — I = 4.08, and a 15 ^m to I ratio 
Si5nm/ Si = 180; the second largest such ratio is 21 for A2218#289. A2219#5 therefore appears 
to be an unusual object, very red in the optical and apparently dominated in luminosity by its 
mid- or far-infrared emission. Lemonon et al (1998) and Altieri et al (1998b) have recently found 
similar objects in ISOCAM images of Abell 2390. 

Referrring to Figure 6, A2219#5 is well-resolved in the optical along the major axis and barely 
resolved along the minor axis. A 2-dimensional gaussian fit gives a FWHM size of 2.0" x 1.0" at 

position angle 103° (major axis). The image seeing is ~ 0.7", so the deconvolved minor axis width 
would be ~ 0.7". The direction to the center of the cluster from the location of A2219#5 lies at 
PA = 210° (see Figure 6), so the major axis is within 17° of being perpendicular to this direction, 
suggesting a lensed arclet. 
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3.3. IR/B Ratios 

Table 1 lists the 15 //m-to-I ratios for all detected galaxies. The three galaxies at the highest 
measured redshifts have the highest ratios, as would be expected for a survey with fewer detections 
in the infrared than the visible. For A2218#289 and A2219#5 at redshifts near 1, the observed 
15 /xm/I ratio maps closely to the ratio of 7 //m/B in the rest frame. For these bands, the observed 
flux density ratios of 21 and 176 translate to a luminosity ratio (in ufi,) of ~ 1 and ~ 10. 

Helou et al (1999) show that for star forming galaxies the ISO 7 /xm band (5 — 9 /xm) is 
dominated by the Aromatic features, and that the luminosity in this band accounts for 2 — 8% 
of the total infrared luminosity between 3 and 1000 /im. The lower luminosity fraction is 
characteristic of intense star bursts, while the higher fractions occur in quiescent or mildly active 
star forming galaxies. At the extreme luminosity and excitation end however, these fractions 
can drop even lower, as in Arp 220, which emits only about 0.3% of its luminosity in the 7 /xm 
band (E. Sturm, private communication; Genzel et al. 1998). Using the higher 7 /xm-to-total-IR 
fraction of 8% leads to IR/B ratios of 13 and 120 for A2218#289 and A2219#5, which are much 
higher than expected in quiescent galaxies. One is therefore led by self-consistency arguments 
to using a smaller fraction applicable to active galaxies. For a 7 //m-to-total-IR fraction of 2%, 
one finds the total IR/B ratios are ~ 40 and ~ 400 for A2218#289 and A2219#5 (extrapolating 
the optical spectrum to obtain Brest), and the total observed IR luminosities about 1.5 x 10^^ 
and 9.5 x W^^Lq. While these arc substantial luminosities, they do have to be adjusted down 
for gravitational Icnsing amplification, which is on the order of ten for A22187^289 (Casoli et al 
1996), and at most a few for A22197^5 given that it is lenscd by the smooth cluster potential. If 
A2219#5 is at a redshift of 0.163 rather than 1.048 (see §2), its luminosity would be considerably 
lower of course but its IR/B ratio would still be very large. 

To further constrain the IR/B ratios we have analyzed the raw IRAS survey data for the 
best estimate of upper limits in the far-infrared, and find for each of A2218#289 and A2219#5 a 
three-sigma upper limit of about 225 mJy at 100 fj,m, in quiet sky with very little cirrus noise. 
This implies that the 7 ^um-to-total-IR fraction must be greater than 1.6% for A2218#289, and 
greater than 1% for A2219#5, or they would have been detected by IRAS. These limit fractions 
add confidence that the 2% adopted above is a reasonably good estimate, since it is constrained 
from both above and below. We therefore conclude that both of the z = 1 objects are likely to be 
dusty galaxies similar to local ULIRGs, with luminosities up to a few x 10^^ Lq. 

While IR/B ~ 40 (A2218#289) is not unusual for an active galaxy, a ratio of 400 (A2219#5) 
is extreme, in that it requires a high column density of gas and dust with the dust surrounding 
the active regions without significant leaks or clumping. Differential lensing amplification between 
wavelengths (Eisenhardt et al 1996) is unlikely to be the cause of the extreme ratio since A2219#5 
is lensed by the diffuse potential of Abcll 2219, which is too smooth to produce such an effect. 
Similar high ratios have been reported by Dey et al (1999) for the z=1.44 ERO (extremely red 
object) HRIO (IR/B ^ 300), and by Yun and Scoville (1998) for IRAS F15307+3252 (IR/B 
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sri 650). 

A2219#5 has a number of other similarities to HRIO and other members of the ERO class, 
although its colors are probably not extreme enough to warrant inclusion in that class (we estimate 
an R— K color of 4.5, interpolating between our measured fluxes; EROs have R— K> 6). The 
EROs are a population of infrared bright, extremely red objects that have been discovered in 
near-infrared imaging surveys. Most EROs are sufficiently faint optically that their redshifts and 
spectral properties are unknown. The exception is HRIO, for which a redshift of 1.44 has been 
derived based on [Oil] and Ha emission lines (Graham & Dey 1996; Dey et al 1999). This object 
is thought to be a high-redshift counterpart of the dusty, ultraluminous infrared galaxies found in 
the local Universe (Dey et al 1999). The optical/IR spectral energy distributions and luminosities 
of A22197^5 (assuming z = 1.048) and HRIO are similar, both are extended at the 1 — 2" level, and 
both show an emission line at long optical wavelengths, firmly identified as [Oil] AA3726, 3729A 
in the case of HRIO and tentatively identified as the same line in A22197^5. If EROs are like 
A221977t5 in their infrared properties, the mid-infrared ISO detection of A2219^5 indicates that 
the sharply rising near-infrared spectrum of objects like HRIO continues out to at least 15 ^m. 
This is consistent with the detection of HRIO at submillimeter wavelengths by Dey et al (1999). 

We conclude that the galaxies in the cluster backgrounds are luminous, high IR/B objects. 
As is the case for most ULIRGs, the primary driving power source, whether QSO or starburst, is 
unclear. 

4. Summary 

Nine galaxies have been detected in ISOCAM 15 /xm images of the lensing clusters Abell 2218 
and Abell 2219. At least three of these are luminous, high IR/B lensed background galaxies at 
redshifts of 0.5 — 1, and one, behind Abell 2219, is an optically faint object heavily dominated by 
its mid-infrared emission. Judging from its high infrared luminosity, its high ratio of infrared to 
optical emission, and its red optical colors, this object appears to be active, with the source or 
sources of activity (starburst and/or AGN) heavily obscured by dust. 

We thank Ian Small for kindly providing optical images and photometry of the clusters, 
Jean-Paul Kneib for the HST image, and Tom Broadhurst and Brenda Frye for obtaining the Keck 
spectrum of A2219#5. 
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FIGURE CAPTIONS 

Fig. 1. — ISOCAM 15 image of Abell 2218. Pixel size is 2", and gray scale range is —40 to +40 
//Jy/pixel. Note coordinate rotation of 40° clockwise (see tick marks). Horizontal striping near the 
bottom of the image is an artifact. 

Fig. 2. — Same as Figure 1, for Abell 2219. Coordinate rotation 44° clockwise. 

Fig. 3. — Contour overlay of ISOCAM image on I-band image from the Palomar 5-m telescope for 
Abell 2218. Dashed contours are negative. Contour levels: (— 3, — 2, — 1, 1, 2, 3, 4)x9 /xJy per 2" 
pixel in original image. 

Fig. 4.— Same as Figure 3, for AbeU 2219. Contour levels: (-3, -2, -1, 1, 2, 3, 4, 5, 6) x7 |LtJy per 
2" pixel in original image. 

Fig. 5.— Contour overlay of ISOCAM image on R-band (F702W) HST image for Abell 2218. 
Contours are the same as in Figure 3. For A2218#289 the centroids of the optical and infrared 
images differ by < 2"; the northwest extending lower contours for this object are caused by the 
image artifact noted in Figure 1. HST image courtesy J.-P. Kneib (see Kneib et al 1996) 

Fig. 6. — Close-up of A2219#5 from the 5-m I-band image (Figure 4). Direction to cluster center 
is almost perpendicular to galaxy major axis, suggesting that this object may be a lensed arclet. 
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Table 1. Galaxies Detected at 15 fim 



Cluster 


ID" 


RA 


Dec 








Sv 


Sb 


Su 


B-I 








(2000) 


(2000) 




(MJy) 


(MJy) 


(MJy) 


(MJy) 


(pJy) 


(mag) 




A2218 


#395 


16 35 48.8 


66 13 02 


0.103 


1100 ± 200 


223 


99 


73 


27 


1.90 


5 


A2218 


#317 


16 35 53.3 


66 12 58 


0.474 


670 ± 200 


45 


11 


6 


2 


2.96 


15 


A2218 


#289 


16 35 54.8 


66 11 52 


1.034 


850 ± 200 


41 


14 


8 


4 


2.52 


21 


A2218 


#275 


16 35 55.4 


66 14 15 




570 ± 200 


66 


23 


16 


10 


2.20 


9 


A2219 


#1 


16 40 13.0 


46 43 07 




890 ± 110 


125 




37 


17 


2.02 


7 


A2219 


#2 


16 40 13.9 


46 41 51 




920 ± 110 


504 




83 


22 


2.65 


2 


A2219 


#3 


16 40 21.5 


46 40 49 




1420 ± 110 


201 




40 


14 


2.45 


7 


A2219 


#4 


16 40 22.5 


46 43 12 




1100 ± 110 


86 




12 


5 


2.79 


13 


A2219 


#5 


16 40 23.0 


46 44 03 


1.048 


530 ±110 


3 




0.13 


< 0.1 


4.08 


176 



^Using numbering system of Le Borgne, Pello, & Sanahuja (1992) for Abell 2218, or RA order for Abell 2219. 

''Redshifts for objects in Aboil 2218 from Ebbcls ct al (1997); redshift for A2219#5 is tentative, based on an emission line 
identified as [Oil] AA3726, 3729A (sec §2 of this paper). 

"^Integrated flux densities of galaxies. Estimated errors « 5% except for A2219#5, where the errors are 8% at I and 15% at 
B. For A2218#289, the photometry covers the main galaxy, and does not include the thin "tail" extending to the northeast 
(sec Figure 5). The /-band photrnctry is Kron-Cousins /, with effective wavelength SOOOA, and an assumed zero point of 2380 
Janskys. The other bands are standard Johnson U, B, V . Photometry courtesy Ian Small. 
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